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ABSTRACT 

Globular clusters (GCs) effectively produce dynamically- formed low-mass X-ray binaries (LMXBs). 
Observers detect ~ 100 times more LMXBs per stellar mass in GCs compared to stars in the fields 
of galaxies. It has also been observationally established that metal-rich GCs are about 3 times more 
likely to contain an X-ray source than their metal-poor counterparts. Recent observations have shown 
that this ratio holds in extragalactic GCs for all bright X-ray sources with Lx between 2 x 10 37 and 
5 x 10 38 ergs _1 . In this Letter, we propose that the observed metallicity dependence of LMXBs in 
extragalactic GCs can be explained by the differences in the number densities and average masses 
of red giants in populations of different metallicities. Red giants serve as seeds for the dynamical 
production of bright LMXBs via two channels - binary exchanges and physical collisions - and the 
increase of the number densities and masses of red giants boost LMXB production, leading to the 
observed difference. We also discuss a possible effect of the age difference in stellar populations of 
different metallicities. 

Subject headings: binaries: close — X-rays: binaries — Galaxies: star clusters: general - globular 
clusters: general 



1. INTRODUCTION 

The preference of bright low-mass X-ray binaries 
(LMXBs), with Lx > 10 36 erg s _1 , to reside in metal-rich 
globular c lusters (GCs) was first noted for our Galaxy 
and M31 (|Grindlavlll99l iBellazzini et al.lll995[ ). Extra- 
galactic observations with a much larger sample of GCs 
and observed bright LMXBs, but also at a higher cut-off 
Lx > 10 37 ergs _1 , have confirmed this ten- 



dency (e.g., 


Kundu et al. 


2002: Maccarone et al. 


2004: 


Jordan et al. 
Sivakoff et a 


2004: Kim et all 120061 Kundu et al. 


2007. 


. 2007). Similar, albeit weaker statistically 



trends have also been seen i n the GCs of M31 for LMXBs 
with Lx > 10 36 ergs~ 1 (|Trudolvubov fc Priedhorskvl 
120041: iPeacock et "ail 120101 ) . As the formation rate of 
bright LMXBs in GCs exceeds that of the field by 
a factor of about 100, it is commonly accepted that 
bright LMXBs in GCs are formed d ynamically (the idea 
was first proposed by I Clarkl 119751 for an overview see 
iVerbunt fc Lewinl I2006TK The origin of the preference 
of bright LMXBs for metal-rich clusters could lie in pre- 
dynamical enhancements of the constituents that eventu- 
ally form into LMXBs, enhancements during dynamical 
LMXB formation, or both. A recent set of deep observa- 
tions of extragalactic LMXBs has strengthened the state- 
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ment that metal-rich cl usters are ~ 3 ti mes more likely 
to have bright LMXBs (|Kim et al.ll2012t) . clarifying that 
this ratio holds across the range of X-ray luminosities 
from 2 x 10 37 ergs _1 to 5 x 10 38 ergs _1 . We show here 
that these latter observations provides the first statisti- 
cally significant constraint on a type of a donor in the 
observed extragalactic bright GC-LMXB population. 

2. DYNAMICAL FORMATION OF LMXBS IN GCS 

There are three sub-populations of LMXBs, depend- 
ing on the evolutionary state of the donor star, that 
can potentially contribute to the X-ray luminosity func- 
tion (XLF) of a dynamically formed X-ray binary (XRB) 
population in a GC: LMXBs with main sequence (MS) 
donors, LMXBs with subgiant or giant donors (hereafter, 
we refer to both giant and sub-giant low-mass stars as red 
giants, RGs), and ultra-compact XRBs (UCXBs) with 
degenerate, white dwarf (WD) donors. 

The vast majority of XRBs in GCs should be ac- 
creting onto neutron stars (NSs) (formation chan- 
nels of XRBs with a black hole (BH) accretor differ 
from those with a NS accretor, see Kalogcr a et al.l 
120041 : Hvanova et al.ll2010l ). which predominantly formed 
through either el ectron-capture supe rnovae or accretion- 
induced collapse (jlvanova et al.ll2008l ) ; the low escape ve- 
locity in GCs leads to the loss of almost all NSs formed 
via normal core-collapse with a large native kick. This 
results in a lighter population of N Ss in GCs, with initi al 
masses of M NS ~ 1.28 ± 0.06 M (jTimmes et al.|[l996l) . 

LMXBs in dense stellar systems are most efficiently 
formed from binar ies with a NS and a MS donor 
(jlvanova et al.l 120081 ) . The strongest formation channels 
among MS donors are binary exchange interactions with 
a NS or the creation of a NS-MS binary following the 
accretion induced collapse of a WD in an already dy- 
namically formed WD-MS binary. Tidal capture is sig- 
nificantly less efficient. Although the formation rate of 
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NS-MS systems is estimated to b e the highest amon g 
the three sub-populations of donors ([Ivanova et al.ll2Q08l) . 
the mass-transfer (MT) rate that these systems with low- 
mass MS donors can drive is very low. As a consequence, 
most NS-MS LMXBs are trans ient with low duty cycles 
and low outburst lumi nosities (jFragos et al.1 120081 120091 : 
iRevnivtsev et al.ll2011|) . As such, they are not expected 
to have any contribution to the observed XLF at the lu- 
minosity range of interest here (Lx > 2 x 10 37 ergs~ 1 ). 
While persistent NS-MS LMXBs can be formed in only 
metal-rich GCs (for old populations systems) and still 
play a significant role for the metal licity dependence in 
our Galactic LMXBs ([ivanov a 2006), their persistent lu- 
minosities are also below the limiting luminosity to which 
extragalactic LMXBs are usually observed. Note that 
the maximum luminosity of 4U 1746-37, the only bright 
LMXB in Galactic GCs that most likely has a MS donor, 
is also below this threshold (jSidoli et al.ll200lh . 

RGs play a significant role in the dynamical forma- 
tion of NS LMXBs, influencing two important formation 
channels. First, RGs can physically collide with NSs 
and form compact binaries that consist of a WD and 
a NS; such binary s ystems may later become UCXBs 
(jlvanova e t al. 2005). Second, following the formation 
(via an exchange encounter) of an eccentric NS-MS bi- 
nary, the MS donor may overfill its Roche lobe as it 
evolves to become a subgiant or giant. In both cases, 
these binaries will appear as persistent LMXBs with 
Lx > 2 x 10 37 ergs~ 1 for a fr action of their MT evo- 
lution (jFragos et al.l[200l I2009T ). 

Dynamical formation of bright LMXBs is determined 
by the probabilities of two events: first, some seed binary 
with a NS has to be formed in a dynamical encounter (see 
§3), and second, this binary must start MT that is fast 
enough to appear as a bright LMXB (see §4). 

3. ENHANCEMENT OF THE DYNAMICAL ENCOUNTERS 

RATE 

The total number of physical collisions (PCs) between 
RGs and NSs is A PC = n RG n NS VRG,NSVoo, where n RG 
and tins are the number densities of RGs and NSs accord- 
ingly, Woo is the relative velocity at infinity and <trg,ns 
is the cross-section of an encounter: 

a / 2G(M RG + M m ) \ 

O-RG.NS = TITp H 5 . (1) 

Here rp is the initial closest approach during an en- 
counter. The number of binary exchange encounters, 
Abe, is calculated similarly, replacing tirg with the num- 
ber density of seed binaries and setting rp » 2a, where 
a is the orbital separation in a pre-exchange binary con- 
taining a RG or a MS star that will evolve into a RG. 

Assuming all stars in a GC are coeval, all RGs in a GC 
have about the same mass, which only varies by a few 
percent compared to the MS turn-off mass of the spe- 
cific GC. This mass range is metallicity dependent since 
the life-times of metal-rich stars are longer than those of 
metal-poor stars, both for their MS life-time, tms, and 
their life-time as a RG, trq- Figure Q] shows the ages of 
a star when it exhausts the hydrogen at its center (end 
of MS) and when it reaches a radius of 30 R Q as a func- 
tion of its zero age MS mass, for two metallicity values 
(Z=0.01 and Z=0.0002), typical of red (metal-rich) and 
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Fig. 1. — Age of star when it exhausts the hydrogen at its center 
(solid lines) and when it reaches a radius of 30 Rq (dotted lines) as 
a function of its zero age main-sequence mass, for two metallicity 
values, Z=0.01 (black lines) and Z=0.0002 (grey lines), typical of 
red and blue GC respectively. 

blue (metal-poor) GCs respectively. The stellar evolu- 
tion tracks used to create Figure [J were calculated u sing 
the MESA stellar evolution code ([Paxton et al.ll20lii ). 

Although the choice of 30 R© in Figured] is not strictly 
rigid, it is justified by two reasons. First, it is approxi- 
mately the maximum radius that a RG can have as part 
of a binary system in a GC environment. The character- 
istic collision time for a binary is 

t co11 » O.SGyr Vl ° (2) 

n5aioo(Mbi n + {M}) 

Here n 5 = n/(10 5 pc~~ 3 ), where n is the stellar num- 
ber density, Mbi n is the total binary mass in Mq, (M) 
is the mass of an average single star in Mq, v\q — 

Voo/ (10 km s~ ) and aioo = a/(100 Rq) where a is the bi- 
nary separation. In a typical dense cluster, a binary with 
a ~ 100 Rq (i.e., about the size of the Roche lobe radius 
for a RG of ~ 30i?©) will be perturbed by a collision dur- 
ing its trg • This likely leads to an increased eccentricity 
that starts MT before the RG reaches ~ 30i? Q . Larger 
orbital separations would correspond to very wide (soft) 
binaries that would be quickly disrupted via strong dy- 
namical interactions. Second, larger RGs are unlikely to 
form a compact binar y after a PC (jlvanova etafl [20051 : 
lLombardi et al.l 120061 ) . Since any RG evolution from 
-Rrg ~ 3Oi?0 through the tip of a RG branch takes less 
than 1% of trg, the exact upper limit of Rq beyond 
3Oi?0 that can make an LMXB/UCXB does not signifi- 
cantly change trq. 

Figure Q] shows that the mass range at which a giant 
star can exist in a metal-rich GC is significantly larger 
than the one corresponding to a metal-poor GC of the 
same age. For a typical age of a Milky Way GC, 11 Gy, 
a red cluster has RGs with masses MRG,red — 0.946 — 
0.972 M , while a metal-poor GC had R Cs with masses 
M RG ,biuo ^ 0.826-0.838 M . Assuming a[Kroupa| ([20011) 
initial mass function, this difference in mass ranges of 
RGs translates to a number density of giant stars tirg 
that is ~ 60% higher in red clusters compared to blue 
clusters. If extragal actic red GC have a younger stellar 
population (see, e.g. lWoodlev et aLll2010l ). the difference 
in number density can be increased even further: e.g., 
RG in a metal-rich GCs of ~ 8 Gyr will be about twice 
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more abundant than in a ~ 11 Gyr metal-poor GC0 The 
total number of dynamical encounters leading to LMXBs 
formation in red GCs is ~ 1.6-2 higher compared to blue 
GCs due to the increase in tirg- 

The average mass of the RGs in a cluster also pro- 
duces a secondary enhancement through cross section of 
encounters, ctrg.ns- Fot PCs, the second term in the 
brackets of Eq. (1) is much larger than one, leading to 
the mass-dependence increase in the cross-sections of en- 
counters in red GCs compared to blue GCs by a factor of 
(M" RG ,rod + M N s)/(M R G,blue + M NS ) . For a typical GC 
age of 11 Gyr, the average mass of a giant star in a metal- 
rich GC is higher (AfaG.red ~ 0.96 M ) compared to a 
metal-poor cluster (Mrcmuc ~ 0.83 M Q ). In a younger 
stellar population, ~ 8 Gyr, MRG.red ~ 1.05 M©. The 
increased cross section provides a mild increase in Npc 
of 5 — 15%. For binary-exchange encounters, this mass 
d ependence is p r esent, but weaker. 

Ilvanova et al.l (|2008l ) found a factor of 3 difference in 
the number of UCXBs formed via PCs and present at 11 
Gyr in GCs that are dynamical twins but have different 
metallicity (see their Table 7, where the ratio is 2.8). In 
that study, which did not separate LMXBs as bright as 
extragalactic LMXBs from quiescent LMXBs that are de- 
tectable only in Milky Way GCs, the importance of this 
channel was not emphasized as it produces a relatively 
small fraction of all LMXBs. However, this channel is 
significant if one only considers bright LMXBs for com- 
parison to extragalactic studies. The reason behind the 
factor of 3 difference was not identified, but re-analysis 
of the simulations showed that this may be due to the 
different number of giants in the stellar populations that 
are present in coeval GCs. The number of RGs at 11 
Gyr in the population of a whole cluster of "standard" 
(metal-rich) model, is 1.85 ±0.05 times bigger than that 
of their "metal-poor" model. At the same time, we found 
that the mass-segregation of RGs in the core may play 
an additional role: in the cores of the clusters, the en- 
hancement of RGs was a factor of 2.1 ± 0.1, likely due 
to different turn-off masses. This suggests the number 
density enhancement in the core of a younger red GC 
may be a factor of 2.6, since mass segregation is likely 
underestimated in these si mulations that are n ot fully 
d ynamically self-co nsistent (|Fregeau et al.l 12009ft . 

iKim et all (|2012l ) also showed that red GCs may con- 
tain more high-luminosity LMXBs (those with Lx > 
5 x 10 38 ergs _1 and thus most likely to have BH accre- 
tors) than blue GCs, with a ratio of 2.5j^ - J. The for- 
mation of BH XRBs differs from NS XRBs '- BH XRBs 
are likely formed via sequences of multiple encounters 
occuring with seed BH-WD binarie s that were formed i n 
an initiating dynamical encounter (|Ivanova "eTan[20T0ft . 
Since PCs with RGs provide at least half of all these 
seed BH-WD binaries, the above metallicity-dependent 
enrichment of RGs will also lead to the production of 
more BH-WD XRBs in red GCs; an exact prediction of 
this overproduction is beyond the scope of current dy- 
namical codes. 

9 The determination of an exact GC age is not possible 
due to a number of obser vational uncertainties (for a review, 
sec Brodic & Strader 2006) and theoretical uncertainties (e.g., 
the presence of bina ries also changes theoretical ages estimate; 
IFan fc de GrTii[20T2T ) . 



4. ENHANCEMENT OF APPERANCE PROBABILITY 
4.1. UCXBs 

The formation of an UCXB via PC between a NS and a 
RG depends not only on the probability of such a collision 
to occur, but also on the capacity of that formed binary 
to shrink sufficiently via gravitational wave radiation to 
start MT. The latter, for each set of companion masses 
and post-collisional eccentricity a nd separation , can be 
determined using equations from iPeterd (|1964f ). Here, 
we denote the maximum post-collisional separation in a 
post-collisional binary as a function of eccentricity that 
is capable of making an UCXB as oucxB(e): only PCs 
where ape < aucxB(epc) lead to UCXB formation. 

The post-collisional binary parameters ape and epc 
are functions of the stellar masses involved in the PC, 
the radius of the RG, -Rrg, the initial closest approach 
during the PC, r p , and Woo (the latter parameter does not 
influence outcomes in GCs). Hydrodynamical studies of 
PCs in GCs between RGs and NSs for various Rg and 
rp have been performed by lLombardi et ail (|2006f ): the 
parametrization of their results, (apc;epc) = F(rp, Rq), 
was used as the underly ing physics in popula tion studies 
of LMXB f ormations in Ilvanova et al.l (|2008f l. 

However. lLombardi et al.l (|2006h orilv investigated two 
RG masses, 0.8 and 0.9 M Q . Thus, any mass dependence 
was not thoroughly explored, although their Fig. 16 sug- 
gests that for larger rp more massive RGs would form 
somewhat closer binaries. Indeed, energy conservation 
implies that the formation of close binaries will be a 
function of the RG's envelope mass M cnv , as the lat- 
ter determines how much energy has to be spend on its 
ejection. Neglecting the kinetic energy at infinity be- 
fore and after the collision, a simplified estimate leads to 
1/opc oc M env . 

Note that RGs located in blue and red GCs of the same 
age would have different envelope masses for the same 
core mass - e.g. for a core of 0.3M Q , M CDV differs by 50%. 
The resulting formed binaries are expected to be more 
compact in the case of a red cluster. Hence RGs in red 
GCs can form a binary (potentially an UCXB) in encoun- 
ters that have a larger initial closest approach njcxB,red 
than RGs in blue clusters. The ratio of these maxi- 
mum closest approaches leading to an UCXB formation 

IS njCXB,rcd/njCXB,bluc ~ Menv.red/Msnv.blue- H cncc i 

the ratio of encounter cross-sections resulting in UCXB 
formation in red and blue clusters is: 



CTJCXB.red ^ ^env.red f -^RG.rcd + A/nS \ ^ ^ g 
CUCXB.blue -^env.bluc \ Aip,G,blue + A/nS / 

(3) 

Combined with § |3l this predicts that metal-rich GCs 
should contain a factor of > 2.7 more UCXBs than co- 
eval metal-poor GCs, in agreement with the results of 
the simulations. Note that since a PC does not lead to 
immediate UCXB formation, the relevant Af env is bigger 
than that of current RGs. However, in all the cases the 
appropriate Af cnv is bigger for metal-rich GCs compared 
to metal-poor ones. 

We tested the expected dependence on the envelope 
mass by performing several PC simulations between a NS 
and RGs of different masses, 0.8 M Q and 1.1 M Q , but the 
same Rb.g and Af coro , using the 3D SPH code Star/Crash 
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Fig. 2. — MT rate (black lines) as a function of time passed from 
the onset of RLOF for three indicative MT sequences. In all case 
the onset of the RLOF happened when the donor star reached a 
radius of ~ 3 R© , and the MT was assumed to be conservative up to 
the Eddington limit. For comparison, the critical MT below which 
the accretion disks become thermally unstable is plotted with grey 
color for each MT sequence. Only the higher metallicity systems 
undergo periods of persistent MT. 



(|Gaburov et al.ll2010t lLombardi et al.ll2011[ ). RGs' stel- 
lar structures were first calculated usin g the STARS/ev 
code and then relaxed in StarCrash (see Glcbbcc k et al.l 
120081 and references therein). The comparison of PC 
outcomes with the same rp but different M cnv con- 
firmed that post-collisional binary separation is a func- 
tion of M cnv in the expected way - smaller for larger 
Mcnv, although the dependence is a bit weaker than a 
strictly linear proportionality (we obtain ratios of post- 
collisional orbital separations being ~ 10% different from 

M e nv,red/ -M enVj bluc) • 

4.2. LMXBs with RG donors 

Given MRCred and MjtG.biue, the mass ratio of the dy- 
namically formed LMXBs at the onset of the Roche-lobe 
overflow (RLOF) will be close to unity (q — M^/M-^s ~ 
0.65 - 0.85). The MT rate that these binaries will drive 
depends strongly on the mass ratio of the binary. Specif- 
ically, the higher the mass ratio, the higher the MT rate 
that the binary will drive. Hence, LMXBs in red GC, 
which have higher mass ratios at the onset of the Roche- 
lobe overflow compared to those in blue GCs, are ex- 
pected to drive overall higher MT rates. 

To more properly study this effect, we used the MESA 
code to calculate three indicative MT sequences between 
a giant donor star and a NS accretor, corresponding to 
a typical LMXB with a giant donor in an HGyr old 
red GC, in an 8Gyr old red GC, and in an 11 Gyr old 
blue GC. Figure [5] shows these calculated MT rates as 
a function of time since the onset of RLOF. The RLOF 
always began when the donor star reached a radius of ~ 
3R0, and the MT was assumed to be conservative up to 
the Eddington limit. For comparison, Figure[2]also shows 
the critical MT belo w which the accreti on disk becomes 
thermally unstable (jDubus et al.lll999h . LMXBs with 



MT rates above the critical MT rate are believed to be 
persistent X-ray sources, while transient LMXBs have 
long-term average MT rates below this critical value. 

Figure [2] shows that the MT rate in the two high- 
metallicity cases reaches peak values of 2.5-5 times 
higher than the low metallicity case, right after the onset 
of the RLO. During this initial phase, LMXBs with gi- 
ant donors become persistent X-ray sources in red GCs, 
while they always remain transient sources in blue GCs. 
As a result, for the same total number of LMXBs with 
giant donors, we expect to be able to detect more bright 
X-ray sources in red GCs than blue GCs. 



5. DISCUSSION 

In this Letter we attempted to explain the observed 
ratio of bright LMXBs in extragalactic GCs of different 
metallicities. We identified that the RG population ca- 
pable of forming close binaries is both more abundant 
(as a fraction among all of the stars) and more massive 
in metal-rich clusters than metal-poor clusters. We pro- 
pose that these properties lead to strong consequences 
on the number of bright LMXBs that can be formed and 
observed in otherwise dynamically similar clusters (total 
GC mass, star's number density and velocity dispersion). 
First, more close binary systems are formed, both via 
physical collision or exchange encounters - the increase 
is a factor of 1.6-2.6 (their fraction among all the stars). 
Secondly, a higher fraction of formed NS-WD will be able 
to start MT as bright UCXBs, as their post-collisional 
configurations are more compact. Finally, LMXBs with 
high-metallicity giant donors drive higher MT rates and 
can appear as persistent systems, while low-metallicity 
LMXBs with giant donors are transient. The combina- 
tion of these factors is capable of producing the observed 
overabundance of bright LMXBs in red GCs compared 
to blue GCs, although a full population synthesis study 
that includes proper physics of both collisions and MT 
with RGs is necessary to quantitatively account for the 
combined effect. 
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